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THE SPECTRUM BETA 
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Beta Lyrae, known eclipsing variable since 
1784, has been the subject large number in- 
vestigations, some which were listed Struve 
his 1957 Russell lecture. The problems posed 
the remarkable spectroscopic behavior this 
system have been discussed Struve [1, 
number occasions, and the first concentrated 
effort solve them was made the Yerkes Observa- 
tory 1940, and simultaneously Kopal 
Harvard. papers were published Volumes 

attempt gather more information Beta 
Lyrae and try answer some the remaining 
puzzles led observational effort Struve and 
Sahade the summer 1955. The new material 
195 plates was obtained Eastman Process 
emulsion with the 32-inch camera the coudé spec- 
trograph the Mount Wilson 100-inch reflector 
which gives dispersion approximately A/mm. 


THE SPECTRUM 


The spectrum Beta Lyrae composite and dis- 
plays absorption well emission features. 
distinguish 

(a) set absorption lines which corresponds 
star spectral type B8. These lines remain reason- 
ably constant intensity and show large variation 
velocity throughout the cycle. 

set peculiar absorption lines, usually re- 
ferred ‘‘the the spectrum,” and es- 


observational material used this investigation was 
obtained Struve and Sahade guest investigators 
the Mount Wilson Observatory. The investigation was partly 
supported grants from the Office Naval Research and the 
National Science Foundation. The authors are indebted the 
Computer Center the University California Berkeley and 
its director, Dr. Henyey, for the use the IBM 701 elec- 
tronic computer. 

John Simon Guggenheim Fellow, February 
1957. 

Numbers brackets indicate references page 32. 


pecially shown and which displays several 
components and must originate expanding shell 
ring which surrounds the entire system. 

(c) Broad and red-displaced absorption lines 
before mid-eclipse, and are usually designated the 
satellite This set can roughly de- 
scribed resembling A2-type star; must origi- 
nate stream gas which moves from the second- 
ary star toward the following hemisphere the 
primary component. 

(d) Violet-displaced absorption lines which appear 
only after mid-eclipse, and are usually designated 
the ‘‘violet satellite These lines display some 
broad feature upon which narrower line 
approximately centrally superimposed, the 
perhaps present only broad absorption. The 
ionization indicated this set lines higher than 
that indicated the red satellite lines. The origin 
the violet satellite lines can traced stream 
from the primary component toward the following 
hemisphere the secondary star. 


TABLE 


USED THE DETERMINATION 
THE RADIAL VELOCITIES 


3679.355 3613.641 4267.15 3831.85 3748.489 
86.833 32.926 3838.39 3736.901 83.347 
91.557 3819.637 3994.996 3933.664 3814.121 
97.154 67.528 4176.164 68.470 3906.037 
11.973 88.646 4390.585 3685.192 4173.450 
21.940 3926.530 4481.228 3741.633 78.855 
34.370 64.727 57.684 4233.167 
50.154 4009.270 3853.657 59.291 96.567 
70.632 26.218 56.021 61.320 4303.166 
97.900 4120.857 62.592 3900.546 51.764 
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TABLE 


RADIAL VELOCITIES THE COMPONENT 


[TRANS. AMER. 


PHIL. SOC. 


Cycle Radial Velocities (Km./Sec.) from 
> 
| | | 
April | | | 
May 
une 
une 
une | 
une 
une | | 


—165.7 —171.1 —168.2 —169.2 

—203.3 —197.7 —191.4 —200.1 

—196.4 —213.3 —203.5 

—200.4 —204.4 —200.7 

—190.1 —177.9 —169.8 —189.2 —182.3 

—186.2 —179.4 —196.3 —180.6 

—144.1 —144.0 —143.7 —137.1 

+115.6 +107.9 +111.3 +112.0 

+128.9 +114.4 +129.3 +117.9 
+124.6 +123.6 +123.6 +120.3 

+123.3 +118.3 +121.7 

+130.9 +125.2 +129.9 +129.7 
+132.5 +132.3 +125.3 +129.6 
+135.8 +128.6 +123.0 +128.6 
+165.9 +159.2 +173.6 +164.5 
+166.8 +162.0 +161.4 +165.7 
+168.6 +163.3 +158.3 +164.7 
+175.7 +169.0 +167.1 
+162.6 +171.3 +171.4 +168.2 
+162.9 +157.4 +168.0 +167.1 
+159.9 +155.0 +163.7 
+161.4 +161.4 +165.2 
+124.6 +127.0 +137.9 +120.5 +130.9 
+115.5 +118.9 +125.8 
+114.3 +115.8 +107.4 +123.7 
+105.8 +115.8 +108.1 +118.8 
+111.2 +104.8 +102.0 +117.7 
+110.7 +107.0 +123.2 

—123.7 —120.8 —116.8 —123.7 

—122.7 —126.0 —126.2 —126.3 

—122.8 —125.3 —124.0 —124.9 

—122.5 —123.1 —125.5 —128.5 

—127.0 —131.9 —137.6 —124.6 

—129.0 —124.0 —130.2 —124.8 

—128.4 —132.0 —132.7 —131.6 —123.8 

—130.0 —132.2 —131.1 —129.0 

—132.2 —135.8 —122.9 —130.2 
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(e) Broad emissions, especially and Hel, 
which give the impression being formed broad 
feature upon which stronger, narrower emission 
superimposed. narrower feature stands out very 
clearly during the eclipse. 


There are striking changes with phase and cycle 
the spectrum. Struve [11] has already described 
some them, but the great complexity the spec- 
trum Beta Lyrae makes impossible give 
description which would justice the wealth 
information contained the material. Therefore, 
atlas appended the end this paper containing 
photographic enlargements our spectrograms the 
region effort has been made line 
the strips, referring them the lines the star. 
Mrs. Zebergs responsible for the photographic work 
involved the preparation the atlas. The phases 
indicated the margins are fractions the period 
and were computed means Prager’s formula [12 
modified Saidov [13], namely, 


The number preceding the decimal point stands for 
the cycle which was arbitrarily assumed zero 
for our first spectrogram, taken April 1955, 
U.T., and corresponding phase 0.0235 
(this spectrogram not reproduced the atlas). 
Table gives the correspondence between calendar 
dates and U.T. mid-exposure and the phases. 


THE RADIAL VELOCITIES 
THE SPECTRUM 


All our spectrograms have been 
radial velocity making settings the absorption 
lines, the laboratory wave lengths and identifications 
which are listed table The results are pre- 
sented table and are plotted figures 

The radial velocities give the most consistent 
set values, with extremely small scatter except 
during the eclipse. Furthermore, they seem 


free from any blending effects. Therefore, they were 
adopted for the determination the orbital elements, 
described later section this paper. The 
velocity curve predicted from these orbital elements 
(listed table last column) shown solid line 
figures 

The velocities from the rest the elements, except 
perhaps and II, and even from the mean 
all elements, depart slightly from the velocity curve 
derived from II, and affect mainly the amplitude. 
probable that these departures are the result 
the streams which are present the system. The 
behavior the radial velocities, especially and 
the region of, and immediately after 
positive velocities, resembles the effect the gaseous 
streams Monocerotis [14] and Coronae 
Borealis [15]. 

several phase intervals the velocities from and 
Hel show systematic trend. For instance, im- 
mediately after secondary minimum the lines the 
Balmer series yield progressively more negative values 
toward the lower members the series; between 
primary minimum and first quadrature both the 
and lines give progressively more positive veloci- 
ties toward the red part the spectrum. Even 
shows the effect the streams and immedi- 
ately after primary minimum the sense that the 
redder the wave length the more positive 
velocity, but there are indications similar effects 
outside primary eclipse. 

The departures the mean radial velocities 
and immediately after secondary minimum are 
worthy further consideration. These means fall 
below the velocity curve from and from other 
elements, and suggest motion away from the star, 
perhaps indicating the existence gases streaming 
out from the system through the external Lagrangian 
point which front the primary component 
the system. This interpretation agrees with the 
spectroscopic behavior the relevant phases. 

The departure the velocities from the 
velocity curve probably spurious: such 
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Radial velocities and velocity curve from II. 
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Fic. Radial velocities from and velocity curve from II. 
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Fic. 10. Radial velocities from and III and velocity curve from IT. 
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Cycle Radial Velocities (Km./Sec.) from Cycle Radial Velocities (Km./Sec.) from 
Date and Date and 
1955 Phase 1955 Phase 
1216 .1051 616 .9693 |—170.2; —106.5; —110.5; —8.7 
951 .2520 —121.5; —72.8 |—135.1 913 .9788 —20.2 +4.1 
—52.6 736 .0509 —117.9 —164.9; —114.0 
900 .4149 —99.1; —32.2 —101.3 602 .1232 —119.6 
734 .6532 —86.2 —99.0 752 .2064 —121.6; —5.2 
806 .7323 —138.6; —82.3 428 .2728 —129.5; —76.0; 
906 .7355 —88.1 —141.0; —82.4 
710 .8066 —89.7 —136.4; —89.4 
1010 .8163 —119.7 28.4 
1158 .8210 —146.8 —127.7 33.8 
600 .8802 |—166.2; —127.7 —122.0 739 .2831 —132.9; —85.7; 
1018 .8940 —129.6 —137.0 432.2 
517 .9552 —127.9 —145.3 
.9614 —131.9 —143.9 24.9 
905 .9675 |—175.4; —130.5 —151.0 +30.4 
.0365 —128.2 —160.3 441 .3509 —133.7; —81.8; —68.6 
.0397 —131.4 —168.1 21.3 
940 .0467 —171.4; —83.4 2.7 
1116 .0518 —132.5 —166.1; —84.3 1.2 
504 .1092 —140.2 —170.5; —118.8 +14.8 
647 .1147 —137.9; —107.3 729 .3599 —134.7; —84.4; 
709 .1159 —128.5; —114.3 28. 
1059 .1282 —138.3; —80.1 —113.2 +16.1 
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TABLE 3—Continued 


Cycle Radial Velocities (Km./Sec.) from 
| +31.8 
+15.9 
+23.8 
+10.6 
| +30.3 
.4357 —135.0; —57.9; —130.1; —77.8; +5.8 
} +6.6 
731 .4373 —132.6; —57.9; —79.9 
+21.5 
| | +25.4 | 
835 .4408 —133.5; —57.2; —78.5 
+30.1 
906 .4425 —131.3; —56.5; —67.8 
| +31.6 
| | | +26.0 
+19.7 
+47.9 
11| .5031 —173.5; —135.8 


+26.8 


+20. 
+26. 
+40.0 
+21.2 
—174.0; —134.7; —46.6; —134.2; —34.2 


+6.1 


+5.2 
+6.0 
1200 .5291 —132.5 —133.6; —29.7 


| 

7| 822 |11.9872 | —-160.8; —122.3; —23.9 —125.6; —35.3 

| 443 .9915 | —163.4; —124.0; —21.0 | —125.0; —30.9 

615 .9965 | —160.3; —124.1; —27.4 | —122.8; —40.2 

8 00 {12.0021 | —165.9; —124.2; —12.5 | —124.4; —7.2 
| | 


departure present when the strong line 
alone plotted. 

During the eclipse the velocities from the lines 
the component show the effect rotational dis- 
turbance. This effect not symmetrical before and 
after mid-eclipse: more pronounced after mid- 
eclipse, fact which could explained terms 
inclination the axis rotation the component 
with respect the plane the orbit. evidence 
for changes the velocity the star either during 
one night different cycles has been found. 


THE SHELL SPECTRUM 


The shell spectrum shows multiple components, 
which some phases blend with the stellar lines. 
The radial velocities from these components 


(of the triplet series) and and 3965 
(of the singlet series), that is, from lines arising from 
metastable levels, are listed table and are plotted 
figures ‘Table lists the radial velocities 
from other lines the triplet series, namely, 
and 4472, and from Hy, which are plotted 
figures and 16. 3888 has two components 
with very definite and rather constant velocities, 
namely, about —170 and —130 km./sec. The rest 
the components not show such constant be- 
havior. The two components which have 
referred must arise shell which surrounds 
entire system. The fact that the more negative 
component observed only immediately 
after mid-eclipse suggests that the density this shell 
relatively higher these phases, conceivably be- 
cause the matter which gives rise the violet 
satellites, part which probably escapes from the 
system through the external Lagrangian point 
front the secondary component the system. 
The different behavior the more negative component 
the lines and also suggests that 
the shell expands outward. 

The shell velocities Mount Wilson coudé 
plates taken Struve 1952 and 1953 are not ap- 
preciably different from those obtained 1955. 
Neither are there any very large differences the 
intensities the lines. would seem that the system 
has been losing mass rather constant rate the 
past few years. 


THE SATELLITE LINES 


The radial velocities from the satellite lines are 
listed tables and and are plotted figure 17. 
the plot the violet satellite lines only the veloci- 
ties from the narrow absorptions have been entered. 
The velocities from the red satellites tend alge- 
braically larger approach mid-eclipse, and the 
opposite seems true for the violet satellites. 


THE EMISSION LINES 


Shifts the position the emissions are roughly 
opposite phase with respect the absorption lines 
the star. The relatively narrow, strong 
the emission structures have been measured for 
radial velocity and the results for 3888 and 4472 
and for are listed table and are plotted 
figures 18, and 20. These velocities are phase 
with those expected for the unobserved lines the 
secondary component the but the y-velocity 
about +130 km./sec. larger than that the 
component. The semi-amplitude the order 
km./sec. 


INTERSTELLAR 


The mean velocity the interstellar and lines 
from 191 plates 15.2 0.2 km./sec. 
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TABLE 


RADIAL VELOCITIES FROM SHELL LINES HEI AND 


Cycle Radial Velocities (Km./Sec.) from 
Date and 
1955 ory Phase 
951 .2520 —84.3 
May 5/1000 2.3408 —35.0; +48.6 —166.7; —41.3; +58.1 
—52.7; +15.7; +67.0 —38.6; +23.5 
703 .4086 —28.7; +35.6 
759 .4116 —28.0; +42.5 —28.0; +45.4 
900 .4149 —23.2; +51.0 —25.9; +39.7 
June 521 4.5687 —141.6; —91.0 —142.7; —70.9; +15.1 
.5725 —142.4; —91.8 —123.8; +5.2 
.5794 —135.0; —66.4 —124.0 
.6595 —142.4; —81.6 —110.2; —34.1 
.6624 —146.9; —103.9 —160.4; —111.0; —64.6 
702 .7288 —97.2; —23.8 —139.8; —51.0 
.7355 —136.3; —73.4; —0.2 —137.2; —47.0 
.7386 —77.6 —130.3; —43.6 
710 .8066 —104.0; —48.7; —5.8 
.8099 —114.6; —64.0 —111.2; —59.8; +3.0 
—117.8; —63.2; —9.8 —119.1; —59.2; —11.2 
.8163 —107.6; —40.1 
600 —121.4; +18.1 —115.5; —2.8 
—119.5 —117.9; —36.4; +29.4 
915 .8906 —122.3; +5.5 —118.7; +4.7 
.8923 —126.0; +32.8 —122.9; +1.1 
.8940 —122.2 —126.5; —19.4 
—129.0; +12.2 —125.1; —6.7 
517 —129.3; —71.2; —8.2 
.9586 —135.4; —67.9 —75.7; +0.4 
.9614 —138.2; —42.4; +20.5 —111.6; —36.9; +23.7 
.9642 —146.8; —75.0; +4.8 —112.4; —33.5; +15.8 
.9675 —134.0; —69.8; +13.5 —116.0; —47.0; +9.4 
.9753 —75.0; —9.8 —81.7; +16.3 
456 —138.5; —92.4 —110.3 
834 .0431 —128.3 —107.7 
504 .1092 —170.6 —154,3 
751 .1182 —170.8 —164.4 
956 —146.5 —165.4 
July 403 6.9621 —100.4; —25.2 —79.6; —15.4 
.9652 —105.8; —19.0 —15.4 
—121.5; —28.5; +18.6 —129.2; —27.7 
.9762 —108.8; —27.2 —125.7; —31.9 
.9788 —115.5; —24.8 —122.2; —70.8; —21.5 


Cycle Radial Velocities (Km./Sec.) from 
Date and 
July 4/ 1002™; 6.9814 | —110.8; —22.1 —97.7; —13.8 
10 47 -9838 —114.2; —57.7; —10.1 —127.3; —70.2; —23.7 
1144 .9869 —122.2; —71.3; —29.4 
§| 355 7.0390 | —120.1 —124.2; —59.3; —4.3 
4 49 .0419 —124.0 —122.8 
5 50 .0452 | —125.4; —60.0 —119.4 
6 41 0480 | —121.1 —118.0 
7 36 0509 | —123.4; —81.1 —110.4 
8 37 0542 —116.6 —118.2 
10 27 .0601 | —117.1 —111.4 
6} 352 -1162 | —70.8 —52.6 
440 .1188 —57.5 —54.8 
§ 21 .1210 | —68.3 —50.6 
6 02 .1232 | —63.4 —54,.9 
6 43 1254 —66.8; +34.2 —52.1; +25.5 
| —171.0; —60.4; +39.4 | —167.0; —57.8; +24.0 
8 08 .1300 | —60.5; +34.2 —65.7 
8 49 .1322 —58.2; +31.2 —57.9 
9 30 .1344 | —161.6; —64.8; +37.7 | —46.0 
10 16 .1368 | —171.3; —61.4; +26.2 | —163.1; —55.3 
1 45 .1400 | —73.0 —47.5 
12 02 —72.3 
.1936 | —125.0; —27.6 —120.5; —36.0 
4 34 1958 | —130.4; —32.6 —122.6; —34.5; +52.9 
5 .1980 | —124.2; —35.6 —115.7; —36.1 
5 56 .2002 —126.8; —23.0 —115.8; —28.4 
6 35 .2023 —123.4; —39.2 —120.7; —15.7 
7 16 .2045 | —121.8; —8.0 —121.5; —25.0 
7 52 .2064 —122.7; —42.7 —122.9; —43.3 
8 27 .2083 —120.0; +4.0 —118.1; —43.4 
9 00 .2101 —123.6; +1.0 —119.5; —7.5 
9 35 .2120 | —119.2; —4.3 —114.7; —25.2 
10 11 .2139 | —118.6; —20.0 —122.4; —24.5 
10 52 .2161 —124.0; —17.7 —117.6; —9.8 
11 38 .2186 —122.5; —1.9 —121.1; —16.1 
8| 3 48 .2707 | —107.8; —1.4 —105.3; +4.6 
4 28 .2728 —111.6; +17.3 —101.8; +12.4 
5 02 .2747 —111.6; +17.1 —103.3; +10.2 
5 33 .2763 | —111.8; +10.3 —103.3; +21.4 
6 02 .2779 | —107.4; +23.5 —99.2; "4-23. 5 
6 35 .2797 —107.8; —0.6 —102.0; +7.9 
7 07 .2814 | —108.0; +10.8 —98.6 
7 39 .2831 —106.2; +18.0 —96.5; +0.8 
8 10 .2848 | —111.8; —6.8 —97.9; +33.2 
8 41 .2864 | —107.0; +12.4 ~110.0; +20.4 
9 11 .2880 | —102.1; —2.3 —94,4; <8 
9 44 .2898 | —100.2; +1.5 —100,9; +3.4 
10 16 2915 —102.8: +3.3 —103.7; +3.4 
10 51 .2934 —105.0; +14.2 —98.2; "+-4.0 
22 .2951 —104,9; +18.4 —100.3; +6.1 
11 53 .2967 | —101 3: +21.0 —96.8; "416.0 
3% 3485 | —77.5; “$2 4 —79.5; +14.2 
441 .3509 | —67.6; +9.7 —61,2; +27.6 
5 26 .3530 | —82.7; +26.7 —70.5; +30.3 
5 53 3547 | —74. 5: +16.2 —72.6; +14.1 
6 30 3567 | —69.4; +10.8 —72.7; +10.5 
7 02 .3585 —76.7; +31.4 —74.1; +9.8 
7 29 .3599 —72.6; +22.8 —55 8: +28.8 
3 00 .3616 | —133. 6; —76.1; +23.8 | —70.7; +13.2 
8 28 3631 —139.7; —65.2; +21.4 | —144.7; —75.6; +4.0 
8 57 3646 | —141.9; —65.5; +12.0 | —58.8; +20.1 
8 29 .3663 | —70.7; "430.4 —141.3; —56.0; +30.0 
10 04 3682 | —132. 4; —75.6; +20.8 | —53.3; ‘S347. 2 
10 36 3699 | —74.9; "421. 5 
11 08 3717 | —75. 7: +26.9 —114.0; —55.5; +31.2 
11 42 .3735 | —72.3; —10.4 —55.2; "4-24. 9 
3 48 4254 | —57.5; +28.6 —59.3; +4.1 
4 28 4275 | —85.9; +28.5 —44.6; +27.3 
5 00 4292 | —70. 6: +13.5 —70.6; +13.5 
5 29 .4308 | —50.9' —47.5 
5 58 4324 | —85.1; +16.1 —47.5 
6 30 4341 —65.3; +20.9 —92.6; —13.7 
7 00 4357 | —77.7; +5.8 —53.3 
7 31 .4373 | —46.2; +6.5 —39.2; +12.3 
8 05 .4392 | —75.0; +5.0 —71.7; +11.5 
8 35 4408 | —73.0 —59.0: —0.5 
9 06 4425 | —82.7; —47.1; +4.9 —76.0; —4.8 
9 35 .4440 | —66.9; +30.2 —4.1 
10 10 4459 —68.4; —9.6 —57.1; +4.9 
10 46 4478 —115.6; —67.0 —122.6; —62.7 
11 22 .4500 —34.6 
11} 11 28 | —141.3 
Sept. 7| 3 22 11. 9872 | —117.7 —90.5 
4 43 9915 | —118.4; —17.3 —99.9; —66.1 
6 15 .9965 | —117.2; —50.4 —85.2 
8 00 /|12.0021 —117.4; —67.2 —101.5; —47.8 
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Fic. 13. Radial velocities from the shell (triplet) line \3888. The distinguish the groupings which the different 
components suggest. Two rather constant velocities about —170 and —130 km./sec. 
components not behave and show large differences different cycles. 
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Fic. 14. Radial velocities from the shell (singlet) lines and 3965. The distinguish the 
groupings which the different components suggest. 
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Radial velocities from the shell (triplet) lines and the cases figures and 14, difficult 


attribute the velocities from the different components different phases, definite groupings. 
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Fic. 16. Radial velocities from the shell line Hy. 
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Fic. 18. Radial velocities from the emission “peak” 
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TABLE 
RADIAL VELOCITIES FROM THE RED SATELLITE LINES 


U.T. Cycle Phase 
June 4.9552 +176.6 
.9586 189.2 

712 192.5 

191.8 

197.5 

184.1 

July 6.9621 171.4 
173.5 

172.9 

185.6 

186.0 

193.6 

194.5 

.9965 215.8 

12.0021 +210.3 

| 


THE LEAST-SQUARES SOLUTION 


Following the conventional notation [16], have 
for the radial velocity the component spec- 
troscopic binary heliocentric system reference: 


The equations condition for the least-squares solu- 
tion are derived taking the variations equation 
(1) was first done Lehmann-Filhés [17]. 
Schlesinger [18] rearranged the equations condi- 
tion form most suitable for practical application 
and also provided the necessary tables [19] facili- 
tate the actual computation the orbital elements. 
The standard form developed Schlesinger still 


TABLE 
RADIAL VELOCITIES FROM THE VIOLET SATELLITE LINES 
Velocity (Km./Sec.) 
Date Cycle and 

1955 Phase (P) Narrow Edges Broad 
Absorption Absorption 
.0431 —408 —144 

July 7.0390 276.5 
.0419 266.6 —356 —153 
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often used. However, some difficulty arises 
Schlesinger’s method when the eccentricity the 
orbit, very small. This difficulty obvious, be- 
cause the time periastron passage, and the 
longitude the periastron, are not well-determined 
quantities for orbits small and they lose their 
gators have arbitrarily fixed applying 
Schlesinger’s method orbits very small eccen- 
tricity and have thereby obtained misleadingly small 
values for the mean errors the orbital elements. 
Luyten [20] proposed 1936 the replacement 
the time nodal passage, which remains significant 
even circular orbit, and Sterne formulated 
1941 the equations condition the least-squares 
solution for the case small eccentricity. Therefore, 
two methods are now available for computing the 
orbital elements: one for large and the other for 
small These are the methods which are best 
suited computer whom only desk machine 
available. 

The introduction the electronic digital computers 
high speed changes the picture. 
Schlesinger’s procedure with the use his many 
tables unnecessary, because the machine can compute 
small fraction second all the quantities that 
are involved the solution. The subdivision into 
two cases the least-squares solution according 
the value can also avoided. Consequently, 
shall reformulate the equations condition 
such way that they will provide most profitably the 
basis for writing (on any electronic digital computer) 
subroutine which will determine automatically the 
orbital elements spectroscopic binary after the 
radial velocities have been fed into the machine. 
essence, the present formulation generalization 
Sterne’s method for orbits very small eccentricity 
form which holds true for orbits any eccen- 
tricity. 

the case Beta Lyrae and also most other 
spectroscopic binaries, the period, well deter- 
mined several series observations which are 
separated time intervals which are long compared 
the period. Therefore, the least-squares solu- 
tion, shall omit the variation leaving alto- 
gether five independent orbital elements that can 
freely varied the solution. From equation (1) 


Following Sterne introduce the 
the mean which the mean longitude 
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TABLE 


RADIAL VELOCITIES FROM THE EMISSION 


[TRANS. AMER. PHIL. SOC. 


June 


Velocity (Km./Sec.) from 


Velocity (Km./Sec.) from 


Cycle and Date Cycle and 
Phase (P) 1955 Phase (P) 
June 
0.0235 +103.3 11516™ 70.4 78.8 
.0999 +101.8 +133.8 +162.5 .0538 62.5 99.9 
+137.1 +161.8 +188.8 85.5 +138.4 +137.1 
+101.0 +143.0 +146.2 
2.3408 +132.3 +181.3 +161.7 98.0 +131.4 +143.4 
+135.3 +161.4 +108.1 +136.1 +155.4 
.4116 +148.8 +157.4 86.8 +119.6 +143.9 
.4149 +117.8 +135.4 +128.5 .1265 +104.1 +138.1 +148.9 
4.5687 +146.3 +141.5 .1299 90.6 +139.3 +138.8 
+146.2 +107.9 +119.3 +104.8 +136.6 +143.8 
.8066 95.1 93.3 +102.2 +130.4 +131.2 
.8099 +105.3 98.9 +110.8 +126.4 +145.5 
.8210 +110.4 99.2 92.8 +140.8 +148.8 
88.4 78.8 89.1 .3599 +168.9 +152.4 +150.7 
.8906 89.9 .3616 +136.6 +156.5 
.9614 82.1 83.5 86.4 4373 +128.5 
.0397 66.0 75.6 5031 +106.5 +100.8 +103.2 


April 
906 
915 
905 
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Fic. 19. Radial velocities from the emission \4472. 
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Fic. 20. Radial velocities from the emission ‘‘peak” Hy. 
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zero. Obviously, remains determinate 
the limit vanishing Since the mean anomaly 
given 


where 

(4) 

have 
Ww 
(5) 


Also introduce two new variables and such that 


Equations (5) and (6) transform the set three 
variables, and into equivalent set three 
new variables: and the latter will used 
the least-squares solution, because they are well- 
defined quantities for all orbits. 

differentiating equations (5) and (6), solving 
for 67, and terms and and sub- 
stituting them into equation (2), obtain the final 
equations condition: 


A by A~6K, +- A36é A46n + A 5019 = Ags, (7) 


where 

and 


Equation (7) has been used the programming 
the solution the machine. Once the equation 
condition derived, can perform the least-squares 
solution according the standard procedure. 

have set the program, the IBM electronic 
data processing machine 701, the least-squares 
solution spectroscopic binary orbit subroutine. 
The only input required this routine, besides the 
observed data (i.e., time, and radial velocity, 
(1) the order magnitude and and (2) 
the approximate time which the radial velocity 
curve reaches its maximum. this way can 
omit the determination preliminary orbit, since 
both the order magnitude and and the time 
nodal passage can estimated from the observed 
velocity curve. 

For each observation, say, the the machine 
solves the Kepler equation internally Newton’s 
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method successive approximations, computes the 
products (z, from 6), adds them the 
previous sums for form 
and then proceeds the next observation. takes 
either part of, the entire, input data ac- 
cording the command fed into the machine the 
start the operation and computes the orbital ele- 
ments. other words, can feed into the machine 
the beginning total number observations 
and use only N’) observations for computing 
the orbit, while the end the computer will print the 
for all observations based upon the final 
solution. After each least-squares solution, the 
are compared internally with the mean errors the 
corresponding quantities: Ay, and 
one more the following inequalities 


én < An, 6T < AT, ( ) 


are not satisfied, another least-squares solution 
performed using the orbital elements determined 
the first solution the preliminary elements for the 
second solution. This process repeats itself several 
times until all inequalities (10) are satisfied. For 
each individual least-squares solution, the computer 
prints: (1) the coefficients the normal equations 
An, the auxiliary orbital elements 
the mean error per plate, and (3) the orbital 
elements Ki, and their corresponding 
mean errors. After the solution which satisfies the 
inequalities given (10) derived, the computer 
automatically shifts the computation the 
all observed points and prints them. also calcu- 
lates and prints ephemeris based upon the final 
orbital elements thirty-three equal intervals the 
eccentric anomaly. Finally, kilometers and 
the mass function units the solar mass are calcu- 
lated and printed. All these calculations and print- 
ings proceed continuously; hence takes only few 
minutes perform these calculations even when 
several hundred observations are used derive the 
elements. Usually, the machine iterates the least- 
squares solution about three six times before all 
inequalities (10) are satisfied. The fixed binary- 
point system has been used throughout the program. 

Using the IBM 701 computer, have computed, 
with the subroutine just described, several sets the 
orbital elements derived from the observational data 
taken different epochs since the turn the century. 
The results are given table together with the 
sources the observational data. Because the 
well-known effect rotational distortion during 
primary eclipse, all observational data taken between 
0.85 and 0.15 were excluded computing the 
orbital elements. When these new sets orbital 
elements are compared with those derived previously 


| 
i 
| 
why 
y 
- 
ial 
i 
= 
fa 
a 
“ 
i 


VOL. 49, PT. 1959] 


THE SPECTRUM BETA LYRAE 


TABLE 


ORBITAL LYRAE FROM RADIAL VELOCITIES DERIVED SEVERAL INVESTIGATORS 


Belopolsky [1] Baxandall [2] Rossiter [3] Sherman [4] 
180.7 188.5 2.0 km./sec. 183.9 0.3 km./sec. 183.6 0.7 km./sec. 
0.026+ 0.032 0.020+ 0.013 0.003 0.018+ 0.005 
10.285+ 0.110 days 0.060 days 0.009 days 9.655+ 0.021 days 
2.642 days 1.908 days 8.601+ 0.350 days 0.787 days 
m.e. per +11.9 km./sec. +10.1 km./sec. 4.4 km./sec. 


After primary minimum. 
References: 


Belopolsky, A.: Astrophys. Jour. 1897. 


Baxandall, E.: Ann. Solar Phys. Obs., Cambridge (part 1): 1-24, 1930. 
Rossiter, A.: Pub. Obs. Univ. Michigan 1934. 


Sherman, F.: Astrophys. Jour. 94: 1941. 


from the same data the early investigators, find 
that some are not good agreement. The disagree- 
ment may attributed two main causes: (1) The 
early investigators may not have excluded the ob- 
servational data taken during primary eclipse. (2) 
error may have been introduced the use 
normal points which were necessary during the era 
desk computers. The first cause more serious 
than the second, can see table where 
have shown three sets orbital elements, all derived 
from the recent Mount Wilson observations. 
first set derived from all observations including those 
during principal eclipse, the second set from all obser- 
vations taken between 12.30 days (0.95 and 0.78 
day (0.06 P), and the third set (the best one) from 
observations taken between 11.0 days (0.15 and 
1.99 days (0.85 P). The difference these three 
sets striking. However, for such small eccentric- 
ity, the value has very little physical significance. 

From tables and (last column) appears that 
except for the y-velocity there systematic change 
with time any the orbital elements. However, 


the change which amounts several km./sec. 
too large attributed chance fluctuation. 
not know present whether this due 
difference the wave lengths used different 
investigators, the existence third body 
the system. Moreover, change with time the 
surface activity the star may also produce varia- 
tion the y-velocity. 


DISCUSSION 


The interpretation the light curve and 
spectroscopic variations Beta Lyrae have been 
discussed several authors, most recently Struve 
the Russell Lecture the American Astronomical 
Society for 1957. Certain aspects this interpreta- 
tion encounter difficulties which have not been re- 
moved our observations. may therefore 
permissible present here some tentative considera- 
tions which view their revolutionary character 
cannot yet regarded sufficiently well established 
displace the earlier interpretation. 


TABLE 
ORBITAL ELEMENTS (AND MEAN ERRORS) LYRAE, BERKELEY (1955) 


Phase interval 
0.007+ 0.002 
9.997+ 0.005 days 
1.425+ 0.676 days 
m.e. per plate 
No. plates 192 


days days 


186.1 185.0 0.3 km./sec. 
0.008+ 0.002 0.017+ 0.002 
9.964+ 0.005 days 10.016+ 0.008 days 
8.875+ 0.495 days 4.884+ 0.256 days 
8.50 
1.9 km./sec. 
155 115 


After primary minimum, determined from the formula used for computing the phases. 
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TABLE the two stars but the side the secondary com- 

VALUES THE MASSES FOR DIFFERENT Mass ponent with respect the center gravity the 
system. The presence the emissions all phases, 

even 0.5 and their behavior during eclipse seem 

153.00 rule out the possibility that the emitting gases are 

1.5 53.1 79.6 located between the two stars and suggest that they 

1.0 surround the secondary component and extend 

distance from its center larger than the radius the 

star. the case Serpentis [26] the broad 

0.3 14.4 4.3 emissions the spectrum shift position with the 

0.2 12.2 2.4 same phase relationship the stellar absorption 
lines, and yield the same velocity amplitude the 


The mass B8V, according the mass luminosity relation, 
about 40. 


Table lists set pairs values for the masses 
the two components, which are consistent with 

Kuiper [9] had reached the conclusion that the 
mass-ratio the order 1.5, implying masses 
about 80© for the primary component and 50© for 
the secondary, and luminosity the order 
for the brighter star. However, Kuiper also remarked 
that the observational evidence favors absolute 
magnitude about —4.5 for this component, value 
which also suggested the number Balmer lines 
which are observed the spectrum. 

accept that the absolute magnitude the 
principal component Beta Lyrae about —4.5 and 
that the mass-luminosity relation holds, the mass 
the star would the order The mass- 
function then yields mass about 22© for the 
secondary component. Consequently, are led 
the conclusion that the primary component Beta 
Lyrae less massive than the secondary. 

Since mass 80© for the primary component 
seems too large for star normal giant 
spectrum and since the luminosity appears 
about —4.5 mag., one wonders whether are not 
dealing with system which the brighter com- 
ponent actually the less massive. possibility 
is, fact, supported the early works Belopolsky 
[23, 24] and Curtiss [25]. 1893 Belopolsky 
measured the radial velocities the broad emissions 
and derived velocity curve opposite phase 
that from the absorption lines the component. 
Belopolsky’s results were confirmed several years 
later Curtiss, who derived semi-amplitude 
km./sec. from the velocities the emissions. 
From our material, Struve [11] estimated that the 
semi-amplitude from the broad emissions the 
order 100 km./sec. 

velocity curve opposite phase that the 
component and with smaller amplitude could mean 
that either the material responsible for the emission 
moves with the secondary star located between 


latter. Moreover, the Wolf-Rayet binaries, the 
motion the component seems 
given the emission lines. therefore reasonable 
believe that Beta Lyrae the velocity amplitude 
from the broad emissions might roughly coincide with 
the velocity amplitude the secondary star the 
system. However, must keep mind that the 
matter which forms the streams and that which con- 
stitutes the expanding shell must contribute the 
broad emissions; therefore, perhaps only spectro- 
photometric study the emission features would 
give information about the relative contributions 
from each source. 

The emission which appear superimposed 
upon the broad emissions and are conspicuous during 
eclipse yield semi-amplitude the order 
km./sec. and y-velocity about 130 km./sec. more 
positive than that the component. The mean- 
ing these results not obvious, but examination, 
especially the region and the agreement 
between the measurements 3888, and 
4472, suggest that they might real. The 
difference then might due the fact that the 
emission are parts Cygni-type profiles, 
the violet absorption edges being undetected because 
the relatively faint continuous spectrum from the 
secondary component. Such conclusion would 
imply that this component abnormal object: 
very massive and underluminous star, which does not 
contribute appreciably the light the system and 
displays observable features only the form emis- 
sions. The satellite lines would still interpreted 
arising from the two streams mentioned earlier. 
The fact that there gaseous stream from the 
secondary star towards the primary must mean that 
the lobe the inner contact surface around the second- 
ary filled with matter. Since the spectroscopic be- 
havior the system suggests that the radius the 
secondary component smaller than that the 
primary, the matter filling the space between the 
secondary star and its lobe the first critical equi- 
potential surface must not very opaque. 

are willing accept that the emission 
and/or the broad emissions originate matter which 
moves with the secondary star, must also conclude 
that the mass the secondary star larger than that 
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the primary. The mass-ratio suggested the 
ratio the velocity amplitude the component 
the velocity amplitude from the broad emissions 
implies (table 10) that the mass the primary com- 
ponent must about and that the secondary 
about 17©. The mass-ratio obtained 
velocity amplitude the emission yields 
mass about for the primary and 11.5© for the 
secondary. 

model Beta Lyrae which the mass the 
secondary star larger than that the primary has 
certain advantages. For one thing, would not need 
attribute very large mass the star. 
could rather attribute this component the mass 
which would expect for normal giant. 

According the current ideas stellar evolution, 
the primary component, being giant, would 
have expanded from its original size the main 
sequence and would now overluminous, terms 
the mass-luminosity relation. The existence 
gaseous stream from the star further suggests that 
probably fills its lobe the inner contact surface. 
the other hand, the secondary component would 
more massive than consistent with the mass- 
luminosity relation. Since now more massive 
than the primary, was almost certainly the more 
massive component when the system was formed and, 
therefore, would have evolved faster than the star. 
Therefore, the former would more advanced 
stage evolution than the latter, being now under- 
luminous for its mass. 

One the puzzles Beta Lyrae which has intrigued 
astronomers for many years the slow velocity 
rotation the component synchron- 
ization between axial rotation and orbital revolution 
required the model Beta Lyrae 
velocity rotation greater than 200 km./sec., the 
profiles the absorption lines show that present 
Vrot km./sec. [27, 28, 29]. Table shows 
that allow the mass-ratio less than one, 
the discrepancy the velocity rotation becomes 
less serious [11, 4]. The remaining difference could 
explained terms evolution [30, since the 


TABLE 


MINIMUM RADIUS AND SYNCHRONIZED ROTATIONAL VELOCITY 
THE COMPONENT FOR DIFFERENT 
ASSUMING THAT THE STAR FILLs LOBE THE 
INNER CONTACT SURFACE 


Minimum Synchronized 
rot 
0.4 200 110 km./sec. 
0.3 
0.15 


The radius B8V star the order 3.50. 
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star has expanded from its original size and now 
fills its lobe the corresponding inner contact surface. 

Since cannot certain whether the motion 
the secondary star shown the broad emissions 
the narrow ‘‘peaks’’ superimposed upon them, 
can more than suggest that the mass the 
secondary may considerably larger than that the 
primary, and that the secondary star certainly and 
the primary quite probably depart from the conven- 
tional mass-luminosity relation. This conclusion, 
far the masses are concerned, not new: was 
recently discussed Gaposchkin [32] (though 
somewhat different manner) and was mentioned 
Struve [11, 4]. 

The light curve Beta Lyrae can satisfied 
variety models which the primary component 
either larger smaller than the secondary. 
impose the condition that the primary star 
smaller mass than the secondary and that fills its 
lobe the inner contact surface, the body which 
front during the principal eclipse will larger than 
the star. Since the spectroscopic observations 
indicate that the size the eclipsing star must 
smaller, new photometric solution would have 
consider the case small secondary star surrounded 
very extended, tenuous envelope which essen- 
tially transparent the continuous spectrum the 
star with its absorption lines. 

The observations Beta Lyrae show that the time 
primary minimum falls somewhat later than was 
predicted the expression which have used 
derive the phases. probably safe derive the 
present value for the correction the predicted times 
minima computing the time superior conjunc- 
tion from the orbital elements table (last column). 
The result 0.379 day 0.0293 later than the pre- 
dicted times minima. Allowing for this correction 
and considering the information contained Struve’s 
1941 paper conclude that the red satellites are 
present least the phase interval 0.92 —0.98 
and that the violet satellites are present the phase 
interval This implies that the 
stream responsible for the red satellite lines seen 
projected upon region the star which located 
closer the radius vector joining the two components 
than the region upon which the stream the violet 
satellites seen projected. The receding (red) 
stream apparently originates the neighborhood 
the Lagrangian point the restricted three bodies 
problem. Since the period rotation the primary 
star larger than the orbital period [3], tempting 
locate the origin the stream which produces the 
violet satellite lines near the tip the primary star 
region slightly displaced the direction opposite 
that the rotation. 
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Following are thirty photographic enlargements 
spectrograms the region 

Table the text gives the correspondence between 
the phases and the calendar dates. 
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